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Abstract 



We develop and discuss in technical detail an infrared-finite factorization and 
optimized renormalization scheme for calculating exclusive processes, which 
enables the inclusion of transverse degrees of freedom without entailing sup- 
pression of calculated observables, like form factors. This is achieved by em- 
ploying an analytic, i.e., infrared stable, running strong coupling as{Q'^) which 
removes the Landau singularity at = ^qcb by a minimum power-behaved 
correction. The ensuing contributions to the cusp anomalous dimension - 
related to the Sudakov form factor - and to the quark anomalous dimension 
- which controls evolution - lead to an enhancement at high of the hard 
part of exclusive amplitudes, calculated in perturbative QCD, while simulta- 
neously improving its scaling behavior. The phenomcnological implications 
of this framework are analyzed by applying it to the pion's electromagnetic 
form factor, including the NLO contribution to the hard-scattering amplitude, 
and also to the pion-photon transition at LO. For the pion wave function, an 
improved ansatz of the Brodsky-Huang-Lepage type is employed, which in- 
cludes an effective (constituent-like) r^uark mass, ruq = 0.33 GeV. Predictions 
for both form factors are presented and compared to the experimental data, 
applying Brodsky-Lepagc-Mackcnzie commensurate scale setting. We find 
that the perturbative hart part prevails at momentum transfers above about 
20 GeV ^, while at lower Q^-values the pion form factor is dominated by 
Feynman-type contributions. The theoretical prediction for the 7*7 — ^ tt^ 
form factor indicates that the true pion distribution amplitude may be some- 
what broader than the asymptotic one. 

ll.lO.Hi, 12.38.Bx, 12.38.Cy, 13.40.Hq, 13.40.Gp 
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Typeset using REVT^iX 



1. Introduction 



The theoretical description of the QCD running couphng ctg {Q^) in the low-momentum 
region has attracted much interest in the last few years [||,|],0,||,||,||J^J^ . In particular, the 
possibility of including power corrections into as (Q^), while preserving its renormalization- 
group (RG) invariance, enables the removal of the ghost (Landau) singularity and restores 
its Q^-analyticity. Such power corrections, sub-leading in the ultraviolet (UV) region, cor- 
respond to non-analytical contributions to the /3-function as to make the running coupling 
well-defined in the infrared (IR) regime but, being not confined within the UV regime, they 
are outside the operator product expansion. 

The existence of power corrections in as (Q^), if true, would greatly affect our under- 
standing of non-perturbative QCD effects. For instance, a power correction to ag gives rise 
to a linear term in the inter-quark static potential at short distances ||^. On a more spec- 
ulative level, one may argue that the source of such terms are small-size fluctuations in 
the non-perturbative QCD vacuum, perhaps related to magnetic monopoles in dual QCD or 
nonlocal condensates. Besides, and in practice, a power-behaved contribution at low scales 
can be used to remove the Landau singularity, present in perturbation theory, supplying 
in this way an IR stable, i.e., ghost-free running (effective) coupling that can be extended 
to the timelike region |pi]|,P|,|3|,^J5|JTI[] . As a result, re-summed expressions, which typically 
involve integrations down to scales /x ^ Aqcd, are not affected by the Landau pole and can 
be safely evaluated. 

The aim of the present work is to develop in detail a factorization and renormalization 
scheme, which self-consistently incorporates such a non-perturbative power correction in 
the running coupling, and then use it to assess and explore exclusive processes. We do 
not, however, propose to involve ourselves in the discussion of whether or not such power 
corrections have a fundamental justification within non-perturbative QCD. We consider the 
ambiguity in removing the Landau pole as resembling the ambiguity in adopting a particular 
(non-IR-finite) renormalization scheme in perturbative QCD. However, this scheme depen- 
dence will be minimized by combining our approach with the Brodsky-Lepage-Mackenzie 
(BLM) commensurate scale setting procedure []T2[. Recall in this context that the parameter 
Aqcd has no special meaning in parameterizing the position of the Landau pole and can be 
traded for an interpolating scale n, on the basis of the renormalization-scale freedom (see, 
for instance, ||^). The justification for such an approach will be supplied a posteriori by the 
self-consistent incorporation of higher-order perturbative corrections and by removing the 
IR-sensitivity of perturbatively calculated hadronic observables. 

A key ingredient of our approach is that the modified running coupling will be taken into 
account not only in the factorized short- distance part, i.e., through the fixed-order perturba- 
tion expansion, but also in the re-summed perturbative expression for the exponentiation of 
soft and collinear gluons (Sudakov effects) and in the RG-controlled evolution of the factor- 
ized parts. This means, in particular, that the exponent of the Sudakov suppression factor 
will be generalized to include power corrections, which encode long-distance physics. 

To accomplish these objectives, we adopt as a concrete power-corrected running coupling 
an analytic model for as, recently proposed by Shirkov and Solovtsov [Q, which yields an IR- 
finite running coupling. This model combines Lehmann analyticity with the renormalization 
group to remove the Landau singularity at = Aqcj-,, without employing adjustable pa- 
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rameters, just by modifying the logarithmic behavior of as by a (non-perturbative) minimum 
power correction in the UV regime. 

At the present stage of evidence, it would be, however, premature to exclude other 
parameterizations for the behavior of the running coupling in the infrared, and one could 
introduce further modifications • It is nevertheless worth remarking that in a recent 



work by Geschkenbein and loffe |jT4| on the polarization operator (related to the Adler 



function) the same infrared limit for the effective coupling was obtained as in the Shirkov- 
Solovtsov approach. Furthermore, it was shown in that postulating that the Adler 



function -D(Q^) is given by integrating R^'^^{s) over the physical region s > only, one 
finds a A-parameterization for the strong running coupling in the spacelike region which 
coincides with the pole-free one-loop expression proposed by Shirkov and Solovtsov [Q. 
Hence, the assumption of analyticity of the strong coupling in the complex plane, used 
by Shirkov and Solovtsov, which at first sight might seem arbitrary, is supported by the 
analyticity of a physical quantity. 

Whatever the particular choice of power corrections in the running coupling, it is clear, 
without the Landau pole, IR sensitivity of loop integrations associated with IR renormalons 
(see |ll| JT5|JT6[| , and for a recent review) is entirely absent. We emphasize, however, 
that the two approaches, though they both entail power-like corrections oc (Aqcd/Q)'', 
are logically uncorrelated, as the removal of the Landau pole is a strong-coupling prob- 
lem (tantamount to defining a universal running coupling in the IR region), whereas IR 
renormalons parameterize in a process-dependent way the low-momentum sensitivity in the 
re-summation of large-order contributions of the perturbative series of bubble chains. In 
some sense the two approaches appear to have complementary scopes: employing "forced 
analyticity" of the running coupling attempts to incorporate non-perturbative input in terms 
of a power-correction term in perturbatively calculated entities, like the hard-scattering am- 
plitude and the Sudakov suppression factor, that are in turn related to observables (e.g., 
form factors). The renormalon technique, on the other hand, tries to deduce as much as 
possible about power corrections from (re-summing) perturbation theory. Whether power 
corrections inferred from renormalons can link different processes (universality assumption) 
is an important question currently under investigation P, p!8| , p!9| , pO| , |2T| . 



Continuing our previous exploratory study p2[ (see also p3[ from which the present in- 
vestigation partly derives), we further extend and test our theoretical framework with build- 
in analyticity by including into the calculation of the pion form factor the next-to-leading 
order (NLO) perturbative contribution to the hard-scattering amplitude [^,^,^,^. To 
compute the pion form factors in the region b ~ I/Aqcdj where the hadronic size becomes 
important, the effects of the original distribution of the partons inside the pion have to 
be included. To this end, an ansatz of the Brodsky- Huang- Lepage (BHL) type for the 
(soft) pion wave function is adopted. This ansatz incorporates an effective (i.e., dynamically 
generated) quark mass to ensure that tp^^^ {x, k_|_) has the correct behavior for k_|_ = and 

— > — oo. Predictions for both the pion electromagnetic and the pion-gamma transition 
form factor, employing a pion light-cone wave function without a mass term, were presented 
in Refs. p2| , p3| to which we refer for further details. The influence of the mass term on the 
hard pion form factor is very weak and it primarily affects the soft contribution and the 
pion-photon transition form factor. Moreover, the specific form of the ansatz for modeling 
the intrinsic transverse parton momentum in the pion bound state is insignificant for the im- 
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plementation of the IR-finite scheme, as its effects become relevant only for large transverse 
distances below 1/Aqcd that are outside the scope of the present investigation. 

To minimize the dependence on the renormalization scheme and scale, we obtain our 
results using an optimized renormalization prescription, based on the Brodsky-Lepage- 
Mackenzie (BLM) [|12| commensurate scale setting. The effect of using a commensurate 
renormalization scale in calculating F^^ and is discussed quantitatively. An important 
observation is that the theoretical prediction for the hard (perturbative) contribution to 
the pion's electromagnetic form factor exhibits no IR-sensitivity, in contrast to approaches 
|,|5D|pT|,|77[] with no IR-fixed point in the running coupling. The existence of an IR-fixed 



point in the effective strong coupling is implied by the general success of the dimensional 
counting rules. 

The major virtue of such a theoretical framework, the latter being the object of this 
paper, is that it enables the inclusion of transverse degrees of freedom, primordial (i.e., in- 
trinsic) and those originating from (soft) gluonic radiative corrections P^ , |S^ , without 
entailing suppression of perturbatively calculated observables, viz., the pion form factor, in 
the high- momentum region, where the use of perturbative QCD is justified, and where sup- 
pression is merely the result of an unnecessarily severe IR-regularization. This enhancement 
is (as noted above) due to power-term generated contributions to the anomalous dimen- 
sions of the cusped Wilson line, related to the Sudakov form factor, and such to the quark 
wave function that governs evolution of the factorized exclusive amplitude. These modified 
anomalous dimensions will be treated here to two-loop accuracy. Note that at the same time, 
the artificial rising trend (see, e.g., [^,0) of the magnitude of the hard pion form factor 
at intermediate and low momenta (below about 4 GeV^), which solely originates from the 
presence of the Landau singularity at = Aqq^) the effective coupling, is here entirely 
absent. Therefore, the scaling behavior of the perturbatively calculated pion form factor 
towards lower values of resembles the one computed with a quasi constant coupling. 
Indeed, the scaled pion form factor is in a wide range of momenta almost a straight line, 
as one should expect for the leading-twist contribution (modulo logarithmic evolutional cor- 
rections which for the asymptotic distribution amplitude start at NLO and are negligible). 
Consequently, the magnitude of the hard part at low is considerably reduced and falls 
short to account for the data. In this momentum regime the pion form factor is dominated 
by its soft Feynman-type contribution. 

Although most of our considerations refer to the pion as a case study for the proposed 
IR-finite framework, the reasoning can be extended to describe three-quark systems as well. 
This will be reported elsewhere. 

The outline of the paper is as follows. In the next section we briefly discuss the essential 
features of the IR-finite running coupling. In Sect. |^ we develop and present in detail 
our theoretical scheme. Sect. | extends the method to the NLO contribution to the hard- 
scattering amplitude. An important ingredient in the phenomenological analysis of the form 
factors is the BHL-type ansatz for modeling [x, k^), which includes an effective quark 
mass. In Sect. |^ we discuss the numerical implementation of our scheme revolving around 
the appropriate kinematic cuts to ensure factorization of dynamical regimes on the numerical 
level by appropriately defining the accessible phase space regions of transverse momenta (or 
equivalently transverse distances) for gluon emission in each regime. In Sect. ^ we apply 
these techniques to the electromagnetic and the pion-gamma transition form factor. We 
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also provide arguments for the appropriate choice of the renormahzation scale and link our 
renormalization prescription to BLM optimal, i.e., commensurate, scale setting. We also 
discuss how our scheme compares with others. Finally, in Sect. |^, we summarize our results 
and draw our conclusions. 



2. Model for QCD running coupling 

The key element of the analytic approach of Shirkov and Solovtsov [^] is that it combines 
a dispersion-relation approach, based on local duality, with the renormalization group to 
bridge the regions of small and large momenta, providing universality at low scales. The 
approach is an extension to QCD of a method originally formulated by Redmond for QED 



At the one-loop level, the Landau ghost singularity is removed by a single power correc- 
tion and the IR- finite running coupling reads 



-an(l)^g2) ^ aPertH) ^q2^ ^ -npert{l) ^g2) 



47r 

Jo 



1 A 

+ 



2 



ln(QVA2) A2-Q2 



(1) 



where A = Aqcd is the QCD scale parameter. 

This model has the following interesting properties. It provides a non-perturbative reg- 
ularization at low scales and leads to a universal value of the coupling constant at zero 
momentum a^^\Q'^ = 0) = 47r//?o — 1.396 (for three flavors), defined only by group con- 
stants. No adjustable parameters are involved and no implicit "freezing", i.e., no (color) 
saturation hypothesis of the coupling constant in the infrared is invoked. 

Note that the IR-fixed point (i) does not depend on the scale parameter A - this being a 
consequence of RG invariance - and (ii) extends to the two-loop order, i.e., a^^\Q'^ = 0) = 
al^\Q'^ = 0) = Q;s(Q^ = 0). (In the following the bar is dropped.) Hence, in contrast to 
standard perturbation theory in a minimal subtraction scheme, the IR limit of the coupling 
constant is stable, i.e., does not depend on higher-order corrections and is therefore universal. 
As a result, the running coupling also shows IR stability. This is tightly connected to the 
non-perturbative contribution oc exp(— 4'7r/a/5o), which ensures analytic behavior in the IR 
domain by eliminating the ghost pole at = A'^, and extends to higher loop orders. Besides, 
the stability in the UV domain is not changed relative to the conventional approach and 
therefore UV perturbation theory is preserved. 

At very low-momentum values, say, below 1 GeV, Aqcd in this model deviates from 
that used in minimal subtraction schemes. However, since we are primarily interested in a 
region of momenta which is much larger than this scale, the role of this renormalization- 
scheme dependence is only marginal. In our investigation we use Aq'^jI^d"^'' = 242 MeV which 

corresponds to Aqcd'^^^ = 200 MeV. 

The extension of the model to two-loop level is possible, though the corresponding ex- 
pression is too complicated to be given explicitly [^. An approximated formula - used in 
our analysis - with an inaccuracy less than 0.5% in the region 2.5 A < Q < 3.5 A, and 
practically coinciding with the exact result for larger values of momenta, is provided by M 
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Po 



1 1 A2 



lnfi + |ln(l + |lnf) 



(2) 



where = H - = 9, /3i = 102 - f nf = 64, and Di = 0.035 for = 3. 

With experimental data at relatively low momentum-transfer values for most exclusive 
processes, reliable theoretical predictions based on perturbation theory are difficult to obtain. 
Both the unphysical Landau pole of as and IR instability of the factorized short- distance part 
in the so-called end-point region are affecting such calculations, especially beyond leading 
order (LO). It is precisely for these two reasons that the Shirkov-Solovtsov analytic approach 
to the QCD running coupling can be profitably used for computing amplitudes describing 
exclusive processes |]3^ , |35|j36[| , like hadronic form factors. The improvements are then: (i) 
First and foremost, the non-perturbatively generated power correction modifies the Sudakov 
form factor [p2| , p7tp8| , p9t^ via the cusp (eikonal) anomalous dimension ||41|], and changes 



also the evolution behavior of the soft and hard parts through the modified anomalous di- 
mension of the quark wave function. This additional contribution to the cusp anomalous 
dimension is the source of the observed IR enhancement (at larger values) of hadronic 
observables and helps taking into account non-perturbative corrections (power terms in Q 
and the impact parameter b) in the perturbative domain, thus improving the quality and 
scaling behavior of the (perturbative) form-factor predictions (at low Q'^). We emphasize in 
this context that the ambiguity of the Landau remover is confined within the momentum 
regime below the factorization scale, whereas above that scale the power correction is un- 
ambiguous. Since we are only interested in the computation of the hard contribution in the 
region ^ k^^, ^ rn^, this ambiguity is in fact of minor importance, (ii) Factorization 
is ensured without invoking the additional assumption of "freezing" the coupling strength in 
the IR regime by introducing, for example, an (external, i.e., ad hoc) effective gluon mass in 
order to saturate color forces at large distances, alias, low momenta, (iii) The Sudakov form 
factor does not have to serve as an IR protector against as singularities. Hence, the extra 
constraint of using the maximum between the longitudinal and the transverse scale as argu- 
ment of as, proposed in and used in subsequent works, becomes superfluous, (iv) The 
factorization and renormalization scheme we propose on that basis enables the optimiza- 
tion of the (arbitrary) constants which deflne the factorization and renormalization scales 



3^ , |37| , ^ , |43[| - especially in conjunction with the BLM commensurate-scale procedure [|12 
This becomes particularly important when including higher-order perturbative corrections 
(see. Sect. ^). 

3. Infrared-finite factorization and renormafization scheme 

Application of perturbative QCD is based on factorization, i.e., how a short- distance 
part can be isolated from the large-distance physics related to conflnement. But in order 
that observables calculated with perturbation theory are reliable, one must deal with basic 
problems, like the re-summation of "soft" logarithms, IR sensitivity, and the factorization 
and renormalization scheme dependence of truncated perturbative expansions. 

It is one of the purposes of the present work to give a general and thorough investigation 
of such questions, as they are intimately connected to the behavior of the QCD (effective) 
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coupling at low scales. 

The object of our study is the electromagnetic pion's form factor in the space-like region, 
which can be expressed as the overlap of the corresponding full light-cone wave functions 
between the initial ( "in" ) and final ( "out" ) pion states: p4| , ^ 

n,\i q ^ 

with 



where the sum in Eq. (J^) extends over all Fock states and helicities Aj (with Cq denoting the 
charge of the struck quark), and where 

1 . ^ / '^-Li + (1 ~ a;i)q_Li , struck quark , . 

1 k_Lj — Xiq±i spectators . ^ ' 

We will evaluate expression using only the valence (i.e., lowest particle- number) Fock- 
state wave function, ipqq{x, k_|_), which provides the leading twist-2 contribution, since higher 
light-cone Fock-state wave functions require the exchange of additional hard gluons and are 
therefore relatively suppressed by inverse powers of the momentum transfer Q"^. Further- 
more, a recent study based on light-cone sum rules, shows that, for the asymptotic pion 



distribution amplitude, the twist-4 contribution to the scaled pion form factor (Q^-^tt (<5^) is 
less than 0.05, whereas the twist-6 correction turns out to be negligible. As we shall see in 
Sect. ^, this higher-twist correction amounts to about 25% of the NLO hard contribution, 
calculated in our scheme. This uncertainty in the theoretical prediction is much lower than 
the quality of the currently available experimental data. 

In order to apply a hard-scattering analysis, we dissect the pion wave function into a soft 
and a hard part with respect to a factorization scale /ip, separating the perturbative from 
the non-perturbative regime, and write (in the light-cone gauge A'^ = 0) 

{x, k^) = rf' {x, kx) e - ki) + (x, k^) d (kl - , (6) 

where the wave function ip^ix, k^) is the amplitude for finding a parton in the valence Fock 
state with longitudinal momentum fraction x and transverse momentum k^ (we suppress 
henceforth helicity labels). Then the large (perturbative) k± tail can be extracted from 
the soft wave function via a single-gluon exchange kernel, encoded in the hard scattering 



amplitude Th, so that ||3^ , |45[ 

^^■-'^ (x, k^) = l'dyj dH^Tn (x, k^; y, U) rf'iv, U) • (7) 
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As a result, the pion form factor in Eq. can be expressed in the factorized form 



€1 ® Cft + Coft 



+ 



soft 



+ 



Th ® Tn^Cft 



(8) 



where the symbol ® denotes convolution defined by Eq. (|^) . The first term in this expansion 
is the soft contribution to the form factor (with support in the low-momentum domain) that 
is not computable with perturbative methods. The second term represents the leading-order 
hard contribution due to one-gluon exchange, whereas the last one gives the NLO correction, 
and the ellipsis represents still higher-order terms. We will not attempt to derive the first 
term from non-perturbative QCD, but we shall adopt for simplicity the phenomenological 

including, in particular, an 



(see also P1 



approach proposed by Kroll and coworkers in ||47[ 
effective (constituent-like) quark mass in the soft pion wave function for the reasons already 
mentioned in the introduction and based on arguments to be given below. This leads to 
a significantly stronger fall-off with of the soft contribution to the space-like pion form 
factor, compared to their approach, whereas the hard part remains almost unaffected - as 
should be the case if factorization of hadronic size effects is preserved. Though Jakob and 
Kroll consider in the option of a Gaussian fc^-distribution with rriq 7^ and argue that 
in that case is significantly reduced, they do not present predictions for the form factor 
and do not follow this option any further (see, however, the predictions in Ref. |^). For 
other, more sophisticated, attempts to model the soft contribution to F^^IQ"^), we refer to 
||,|8|g|,|50|. 

We now employ a modified factorization prescription [^pO| , which explicitly retains 
transverse degrees of freedom, and define (see for illustration Fig. Il|) 



^f' (ki < §) exp 




ki < cU'Q 



THard > ki > 



(9) 



with b, the variable conjugate to k_i_, being the transverse distance (impact parameter) 
between the quark and the anti-quark in the pion valence Fock state. The Sudakov-type 
form factor exp (—5) comprises leading and next-to- leading logarithmic corrections, arising 
from soft and collinear gluons, and re-sums all large logarithms in the region where Aq^j-, <C 
ki <^ l^;^;^- The source of these logarithms is due to the incomplete cancellation 
between soft-gluon bremsstrahlung and radiative corrections. It goes without saying that the 
function S includes anomalous-dimension contributions to match the change in the running 
coupling in a commensurate way with the changes of the renormalization scale (see below 
for more details). 

Going over to the transverse (or impact) configuration space (typical in eikonalization 
procedures), the pion form factor reads [p9| a 



1 foo /72u 

dxdy / — Vr {y, b, P'; C,, C2, C) Th (x, 6, Q- C^, C) 

J-00 [Any 

X P^(x,6,P;Ci,C2,C4) , 



(10) 



Note that this expression cannot be directly derived from Eq. (P). 
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Figure 1. Illustration of the factorized pion form factor, exhibiting the different regimes of dynamics. 
The wiggly line denotes the off-shell photon. Gluon exchanges are not explicitly displayed. The region 
of hard-gluon re-scattering (LO and NLO) is contained in the short-distance part, termed Th. The blobs 
e{—S) represent in axial gauge Sudakov-type contributions, whereas non-perturbative effects are absorbed 
into the (universal) pion wave functions and ^/^°"*. 

where the modified pion wave function is defined in terms of matrix elements, viz., 

K(x,6,P,/i) = y'"'"'^''ci2k^e-*'^^-^^K(x,k^,P) 

= / ^e-^"-(0|T(g(0)7+75g(0,^-,b^))|7r(P))^^^^ (11) 

with P+ = = P^', = — (P' — Py, whereas the dependence on the renormahzation 

scale /i on the rhs of Eq. ( |Tl| ) enters through the normalization scale of the current operator 
evaluated on the light cone and the dependence on the effective quark mass has not been dis- 
played explicitly. Note that in the light-cone gauge A'^ = 0, the Schwinger string in Eq. (|1T]) 
reduces to unity. The factorized hard part Th (x, y, b, Q; C3, C4) contains hard-scattering 
quark-gluon subprocesses, including in the gluon propagators power-suppressed corrections 
due to their transverse-momentum dependence. These gluonic corrections become impor- 
tant in the end-point region {x 0) for fixed Q^. Furthermore, current quark masses, being 
much smaller than the resolution scale (set by the invariant mass of the partons) can be 
safely neglected in Th, so that (valence) quarks are treated on- mass shell .i 

A few comments on the scales involved and corresponding dynamical regimes (see 
Eq. (|10D and Fig. ^ are in order: 

• The scale C3/6 serves to separate perturbative from non-perturbative transverse dis- 
tances (lower factorization scale of the effective sub-sector and transverse cutoff). We 
assume that some characteristic scale &nonp — 0<-±y^'^/x{l — x) ~ 0.5 GeV exists, re- 
lated to the typical virtuality (off-shellness) of vacuum quarks. This scale should also 
provide the natural starting point for the evolution of the pion wave function. In the 
following, we match the non-perturbative scale C3/6 with the scale Ci/b, where the 
re-summation of soft gluons in the effective sub-sector starts, i.e., we set Ci = C3. The 



^Furthermore the chiral limit is adopted here, i.e., = 0, since the pion mass is much smaller 
than the typical normalization scale in Eq. (11). 
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Figure 2. Regimes of chromodynamics characterized by scales typical for the quanta involved. 

lowest boundary of the scale Ci/b (IR cutoff) is set by Aqcd, though the results are 
not very sensitive to using a somewhat larger momentum scale, as we shall see later. 

• The re-summation range in the Sudakov form factor is limited from above by the scale 
C2^Q (upper factorization scale of the effective sub-sector and collinear cutoff ).i This 
scale may be thought of as being an UV-cutoff for the effective sub-sector, i.e., for the 
Sudakov form factor, and enables this way a RG-controlled scale dependence governed 
by appropriate anomalous dimensions within this sub-sector of the full theory. 

• Analogously to these factorization scales, characterized by the constants Ci,C2, and 
C3, we have introduced an additional arbitrary constant C4 to define the renormaliza- 
tion scale C^fi^x, y)Q = /xr, which appears in the argument of the running coupling 
a'^^ (choice of renormalization prescription). This constant will play an important role 
in providing the link to the BLM (commensurate) scale-fixing. The running coupling 
plays a dual role: it describes the strength of the interaction at short distances (in 
the fixed-order perturbation theory), and controls via the anomalous dimensions of 
the cusped Wilson (world) line and the quark field, respectively, soft gluon emission 
and RG-evolution of Th and to the renormalization scale. The important point 
here is that the analytically improved running coupling contains a nonperturbative 
contribution, which reflects the nontrivial structure of the QCD physical vacuum. 

The appropriate choice of the unphysical and arbitrary constants Cj will be discussed in our 
numerical analysis in Sect. ^. 

The ambiguities parameterized by the scheme constants Cj emerge from the truncation of 
the perturbative series and would be absent if one was able to derive all-order expressions in 
the coupling constant. In fact, the calculated (pion) form factor depends implicitly on both 
scales: the adopted renormalization scale via ctg, and the particular factorization scheme 
through the anomalous dimensions. Since the latter also depend on as, the factorization- 
scheme and the renormalization-scheme dependences are correlated. On the other hand. 



^Note that the constant C2 here differs in notation by a factor of \/2 relative to that used by 
Collins, Soper, and Sterman in (H), i.e., = V^C2- 
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the physical form factor is independent of such artificial scales and satisfies fi — ^^-^ — - = 0, 
for /i being any internal scale. Obviously, both scheme dependences should be treated 
simultaneously and be minimized in order to improve the self-consistency of the perturbative 
treatment. In order to render the perturbative prediction reliable, the parameters Ci should 
be adjusted in such a way, as to minimize the influence of higher-order corrections, thus 
resolving the scheme ambiguity. However, in the present investigation we are not going 
to explicitly match the fixed-order NLO contributions with the corresponding terms in the 



re-summed expression for the "soft" logarithms. Available calculations p^j25|j26|j27[] of the 
NLO contribution to the hard-scattering amplitude for the pion form factor do not include 
the k_i_ components of the gluon propagators, making such a task difficult for the moment. 
Instead, we are going to show in Sect. |^ that a potential double counting of re-summed and 
NLO contributions is de facto very small and of no real practical importance, especially in 
view of the poor quality of the existing experimental data. In addition, to limit a possible 
double counting as far as possible, we have meticulously restricted the numerical evaluation 
of our analytic expressions to the appropriate kinematical regimes. 

In Eq. (|T0|) , Th is the amplitude for a quark and an anti-quark to scatter via a se- 
ries of hard-gluon exchanges with gluonic transverse momenta (alias inter-quark transverse 
distances) not neglected from the outset. To leading order in the running coupling, one has 

Th (x, y, b, Q; /.r) = SCpaTifiDKo (^bQ) . (12) 

This result is related to the more familiar momentum-space expression 

TH(x,i/,kx,U,g,/iR) = — — — -—2 (13) 

xyQ^ + (k_L + l±) 

via the Fourier transformation 

/oo 
d'^hx Th (a;, y, 6, Q, hr) exp [ib_L ■ (k_L + l±)] , (14) 
-oo 

where use of the symmetry of tp.,^ under x ^ 1 — x = x has been made, and where Cp = 
[N^ — l)/2Nc = 4/3 for SU (3)c. In the limit Q"^ oo (and x fixed) this expression coincides 
with the asymptotic hard scattering in the coUinear approximation, up to suppressed power 
corrections. The latter become important in the end-point region {x 0) at fixed Q^, where 
the actual momentum fiow in the gluon propagator becomes small (typically of the order of 
^qcd) so that gluonic transverse momenta cannot be safely neglected. 
The amplitude 



Vn {x, b,P c::^Q, Ci, C2, fJ.) = exp 



-s {x, b, Q, Ci, C2) - s {x, b, Q, Ci, C2) 



2 r -i,{aTm 

JCi/b fl 



VA^ACi/b) (15) 



describes the distribution of longitudinal momentum fractions of the qq pair, taking into ac- 



count the intrinsic transverse size of the pion state |30] and comprising corrections due to soft 



real and virtual gluons including also evolution from the initial amplitude (x, 6, Ci/b) 
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at point Ci/h to the renormahzation scale (x Q. Let us emphasize at this point that the 
power-behaved term, a^^*^^^, does not change the leading double logarithmic behavior of the 
Sudakov exponent. The main effect of the absence of a Landau pole in the running coupling 
a™ is to make the functions s (x, 6, Q, Ci, C2), s (x, 6, Q, Ci, C2) well-defined (analytic) in 
the IR region and to slow down evolution by extending soft-gluon cancellation down to the 
scale Ci/6 ~ Aqcd, where the full Sudakov form factor acquires a finite value, modulo its 
dependence (see Fig. 0). In addition, as we shall see below, the Sudakov exponent con- 
tains power-behaved corrections in {Ci/hKf'^ and (C2/^QA)^^, starting with p = 1. Such 
contributions are the footprints of soft gluon emission at the kinematic boundaries to the 
non-perturbative QCD regime, characterized by the transversal (or IR) and the longitudinal 
(or coUinear) cutoffs. 

The pion distribution amplitude evaluated at the (low) factorization scale Ci/h is ap- 
proximately given by 

(x, 6, Ci/h, niq) ~ ^ ^' "^'^^ ' ^^^^ 



' c 

Because we retain the intrinsic transverse momenta of the pion bound state, we have to 
make an ansatz for their distribution. In the present work, we follow Brodsky, Huang, and 
Lepage (see, also |l30| , |47|| ) and parameterize the distribution E(x, k_L, mq) in the intrinsic 
transverse momentum k±_ (or equivalently the intrinsic inter-quark transverse distance h) in 
the form of a non-factorizing in the variables x and k± (or x and b) Gaussian function, so 
that 

(x, k^, mq) = M^^^{x, Ci/6)S (x, k^, mq) , (17) 



where 



<l>,(x,Ci/6) = A<l>as(x) = A6x(l-x) (18) 



is the asymptotic distribution amplitude, with A being an appropriate normalization factor, 
and where 

S (x, k_L, mq) = 167r^/9^(y'(x)S (x, k^) S (x, mq) (19) 

with 

t (x, kx) = exp \-j3V9{x) \ , (20) 



and 



E(x,mq)=exp -I3'^m\g{x) , (21) 



models the distribution in the intrinsic transverse momentum in the form of a Gaussian in 
the sense of the BHL ansatz.l 



^The width parameter (5 of the Gaussian distribution should not be confused with the /3-function. 
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Neglecting transverse momenta in Eq. (|T3D (collinear approximation), the only depen- 
dence on k_i_ resides in the wave function. Furthermore, limiting the maximum value of k±, 
these degrees of freedom can be integrated out independently for the initial and final pion 
states to give way to the corresponding pion distribution amplitudes, which depend only 
implicitly on the cutoff momentum /x^ = (Ci/6)^: 



(pT, {x, Ci/h) 



167r3 



{x, k_L, m, 



(22) 



where = 130.7 MeV and Nc = 3. Setting g{x) = 1/ (a;(l — x)) and integrating on both 
sides of this equation over x, supplies us with the constraint 



1 = A dx6x{l — x) exp 
Jo 



x(l — x) 



(23) 



because by virtue of the leptonic decay vr /i^t'/i the rhs of Eq. ( p2D is fixed to ^^ifw 
that 



dx (a;, /i^, = 



(24) 



and 



167r3 



S (x, kj 



(25) 



for any factorization (normalization) scale /i. The contributions from higher Fock states are 
of higher twist and contribute corrections at higher order in 1/Q, i.e., they are power-law 
suppressed. 

Moreover, from 7r° 'j'-f, we have 



1 = dx exp 

Jo 



x{l — x) 



(26) 



while the quark probability Pqq and the average transverse momentum (k^^) are given by 



Pqq = 12AV>^/5^ C dx^i^ - exp 
Jo 



x{l — x) 



(27) 



(k^2) = / dxx^il - xfexp 

Pqq Jo 

By inputting /,r and the value of the quark mass mq, we determine the parameters A, (3"^, 

Pqq, and (k_L^)^^^ (see. Table |1]). 

How large should be the quark mass used in the BHL ansatz? The parameter mq has not 
the meaning of a real mechanical mass for the quark, but reflects the complicated structure 
of the QCD vacuum. Let us make this point more transparent. The QCD Lagrangian 
contains no mass scale in the chiral limit. A mass scale enters perturbatively only through 



2/5^m, 



2 



x(l — x) 



(28) 
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Table 1. Values of parameters entering the pion wave function, using the notations of [30|. The 
values in parentheses refer to the case mq = 0. Here the subscript as on refers to the asymptotic 
distribution amplitude. 



Input parameters 




Determined parameters 


mq = 0.33 GeV 




^ = i- 10.01 (i-6) 


U = 0.1307 GeV 


Pi 


= 0.871 GeV-2 (0.743 GeV'^) 






2)^/2 ^ Q_352 (0.367 GeV) 






Pqq = 0.306 (0.250) 



dimensional transmutation to enable renormalization. Nonperturbative scales derive from 
vacuum fluctuations of some deflnite correlation length in the context of specifled QCD 
vacuum models. For instance, in the chiral quark model derived from the instanton approach 
( ||5^J53[] ), the pivotal nonperturbative scale is the average instanton size p, whose inverse 
deflnes a mass scale of the order of 0.600 GeV. Breaking chiral symmetry spontaneously in 
the instanton vacuum by the delocalization of the fermionic zero modes, the massless quark 
acquires a momentum-dependent mass to become a quasi-particle. The obtained value of 
this effective mass in the quark propagator is M ~ p/i?^ ~ 0.300 — 0.350 GeV, where R 
denotes the separation between the quark and the antiquark. Note that though one deals 
with massive quarks, higher Fock states are not necessarily zero, as the quark propagator 
contains a renormalization factor Z ~ 1 + 0{p'^ /B?) from which one infers that parton 
quarks and the "constituent" quarks (the quasi-particles) of this model are equal at leading 
order in the small parameter p^/R^. Hence, on the basis of the nonperturbative structure 
of the QCD vacuum, we may conclude that the mass scale characterizing the quarks in the 
pion is of the order of the typical constituent quark mass. Then, it seems plausible to set 
the mass scale mq in the BHL ansatz equal to the dynamical mass Mq obtained in such a 
nonperturbative vacuum model, rather than to use a current quark mass of a few MeV. Only 
for a very dilute instanton vacuum one can realize Mq -C p^^, but then the shape of the pion 
distribution amplitude deviates signiflcantly from the asymptotic one becoming very flat. 
Actually, one may considerably simplify the BHL ansatz by setting all scales responsible 

1 /2 

for the intrinsic (transverse) structure of the hadron (the pion), namely (k_L) and mq 
equal to Aqcd which is of the same order of magnitude: 0.200 — 0.350 GeV.Q Indeed, the 
predictions obtained with this simpliflcation are close to those calculated with the values 
given in Table |1|. 

Let us now return to Eq. ([T5|) . The Sudakov form factor Fs (^, 6, Q, Ci, C2), i.e., the 
exponential factor in front of the wave function, will be expressed as the expectation value 
of an open Wilson (world) line along a contour of flnite extent, C, which follows the bent 
quark line in the hard-scattering process from the segment with direction P to that with 
direction P' after being abruptly derailed by the hard interaction which creates a "cusp" 



^The average transverse momentum of the pion was recently determined in [50| on the basis of 
local duality. Values between 0.260 and 0.320 GeV were found, obviously consistent with the actual 
value of Aqcd and those in Table ||. 
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in C, and is to be evaluated within the range of momenta termed "soft", confined within 
the range hmited by Ci/h (IR cutoff) and C2iQ (longitudinal cutoff) (where ^ = x, x, |/).i 
Thus we have 0,|^,|9|,|O|,g 

Fs {W{C)) = ^Pexp [tg j^dz ■ t'' A^z)^! ^^^^ , (29) 

where P stands for path ordering along the integration contour C, and where < ... >a denotes 
functional averaging in the gauge field sector with whatever this may entail (ghosts, gauge 
choice prescription, Dirac determinant, etc.). Having isolated a sub-sector of the full theory 
(cf. Fig. where only gluons with virtualities between Ci/h and C2^Q are active degrees 
of freedom, quark propagation and gluon emission can be described by eikonal techniques, 
using either Feynman diagrams or by employing a world-line casting of QCD which 

reverts the fermion functional integral into a first-quantized, i.e., particle-based path integral 

Then the Sudakov functions, entering Eq. (p!5[), can be expressed in terms of the 
momentum-dependent cusp anomalous dimension of the bent contour to read 

s (e, 6, Q, Ci, C,) = - r,,,p (7, aTif^)) (30) 

2 Jc'i/b fl 

with the anomalous dimension of the cusp given by 

rcu.p(7,<(/i)) =21n(^^^A(ar(/x))+i?(ar(/^)) , 

ppcrt _i_ pnpcrt /n-i "\ 

cusp cusp 5 v"-* / 

7 = ln(C2^Q//i) being the cusp angle, i.e., the emission angle of a soft gluon and the bent 
eikonalized quark line after the external (large) momentum Q has been injected at the cusp 
point by the off-mass-shell photon, and where in the second line of Eq. ( ^T]) the superscripts 
relate to the origin of the corresponding terms in the running coupling (cf. Eqs. (|1]), (^). 
The functions A and B are known at two-loop order: 



7^(ar(/x))+/3(^7)^/C(Ci,ar(/^)) 



and 



TT 2 \ TT / 



B (ar (/i)) = -I [JC (Ci, ar (/i)) + G iC, C^, aTifi))] 



(32) 



^This means that the region of hard interaction of the Wilson line with the off-shell photon is 
factorized out. 
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The first term in Eq. ([3^ ) is universal,! while the second one as well as the contribution 
termed B are scheme dependent. The K-factor in the MS scheme to two-loop order is given 

by pjsajigpg 
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rifTp + Z^oln (Cie^V2) 



(34) 



with Ca = Nc = 3, rif = 3, Tp = 1/2, and 7e being the Euler-Mascheroni constant. 

The quantities /C, Q in Eq. (|33|) are calculable using the non-Abelian extension to QCD 
4^ of the Grammer-Yennie method |^ for QED. Alternatively, one can calculate the cusp 
anomalous dimension employing Wilson (world) lines p8| , |39| , |40| , |5^ .l In this latter approach 
(see, e.g., |0), the IR behavior of the cusped Wilson (world) line is expressed in terms of 
an effective fermion vertex function whose variance with the momentum scale is governed 
by the anomalous dimension of the cusp within the isolated effective sub-sector (see Fig. |]). 
Since this scale dependence is entirely restricted within the low-momentum sector of the full 
theory, IR scales are locally coupled and the soft (Sudakov-type) form factor depends only 
on the cusp angle which varies with the inter-quark transverse distance b ranging between 
Ci/6 and C2^Q. 

The corresponding anomalous dimensions are linked to each other (for a nice discus- 



sion, see ||3^) through the relation 2r 



cusp 



(a: 



7yc(ar(/^)) with T 



cusp 



.a: 



Cf a^^{fi^)/Ti, which shows that ^jjc = A{a^^{fx)). (Note that 7g = — 7/c-) 

The soft amplitude V-j^ (x, b, Ci/b, ^) and the hard-scattering amplitude Th (x, y, b, Q, /i) 
satisfy independent RG equations to account for the dynamical factorization (recall that 
both b and ^ are integration variables) with solutions controlled by the power-term modified 
"evolution time" (see, e.g., and earlier references cited therein): 



r 



In (/iVA') 



/?o In (clI ibk)' 



1 



In- 



{Cr/by 



In 



1/^' - A' 



9i 

62 



A2 



(35) 



from the factorization scale Ci/b to the observation scale yU, with A denoting Aqcd as 
before. The evolution time is directly related to the quark anomalous dimension, viz., 
7q(as'^(/i)) = — ag°(/i^)/7r. One appreciates that the second term in ( p5|) stems from the 
power-generated correction to the running coupling, a^P''''*, and is absent in the conventional 
approach. At moderate values of /i^ this term is "slowing down" the rate of evolution. 



^In works quoted above, the cusp anomalous dimension is identified with the universal term, 
whereas the other (scheme and/or process dependent) terms are considered as additional anomalous 
dimensions. Here this distinction is irrelevant. 

^The derivation of the cusp anomalous dimension in the l/N^ approximation (single-bubble-chain 
approximation) was given in ||2^, Appendix A. 
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The leading contribution to the IR- modified Sudakov functions s b, Q, Ci, C2) (where 
^ = x,x, y, y) is obtained by expanding the functions A and 5 in a power series in a^^ and 
collecting together all large logarithms In (^^^Q) 5 which can be transformed back 

into large logarithms In ^Q^/k^^) in transverse momentum space. Employing equations (|ip 
and d^), the leading contribution results from the expression 

s(i,^Q,C,,C.) ^ \ j2.„ f^) [<i!M^<.. + (CO 

2 Jci/b ^ \ \ ^ ) \_ i: \ TT J 

+<^MBa.(c,c,) + of<^)n. (36) 



7r V TT 



where Eq. is to be used in front of A^^^ , whereas the other two terms are to be evaluated 
with Eq. (0). The specific values of the coefficients A^^\ B^'^^ are 




(37) 

in which the term proportional to A^^'' represents the universal part. As now the power- 
correction term in a^^ gives rise to poly-logarithms, a formal analytic expression for the full 
Sudakov form factor is too complicated for being presented. We only display the universal 
contribution in LLA: 




(38) 



where Q represents the scale C2^Q and the IR matching (factorization) scale /ip varies with 
the inverse transverse distance b, i.e., /ip = Ci/b. Note that the four last terms in this 
equation originate from the non-perturbative power correction (cf. Eq. (|3TD), and that Li2 is 
the dilogarithm (Spence) function which comprises power-behaved corrections of the IR (6A) 
and the longitudinal (Q/A) cutoff scales. In the calculations to follow, Eq. (^) is evaluated 
numerically to NLLA with appropriate kinematic bounds to ensure proper factorization at 

the numerical level. Above (see Eq. (|36D ), we have replaced (a^^^^ by (^a^^^^^^ . Here 
we have an analytization ambiguity. Since nonlinear relations are not preserved by the 
analytization procedure (see, also |ll[), we could have made the square of the running 



coupling, (ai^^^ , analytic as a whole. We plan to report on these interesting conceptual 
issues of analytization in a separate publication. 

Note that, neglecting the power-generated logarithms, we obtain an equation for the 
conventional Sudakov function, which we write as an expansion in inverse powers of the first 
beta-function coefficient jSo to read 
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s{^,b,Q,C,,C2) 



1 

Jo 



(2A«Q + 5«)ln^-2A« [Q-b 



Pi 
Pi 



Q Q-b 



1 + In 26 



1 r. 



In^ (2Q) - In^ (26)1 ' (39) 



where the convenient abbreviations [^] Q = ln and 6 = In ^ have been used. 

This quantity differs from the original result given by Li and Sterman in ^9] , and, though 
it almost coincides numerically with the formula derived by J. Bolz it differs from that 
algebraically. 

All told, the final expression for the electromagnetic pion form factor at leading pertur- 
bative order in Th and next-to-leading logarithmic order in the Sudakov form factor has the 
form 



F^iQ') = - A^Cf fl / dx dy bdbaT'-^^ (^r) $as(x)$as(y) exp 
O JO Jo Jo 



6^ {xx + yy) 



X exp 



^ \xx yy 



where 



Ko iV^Qb) exp [-S {x, y, 6, Q, C,, C2, C,)] , (40) 



S {x, y, 6, Q, Ci, C2, C4) = s {x, 6, Q, Ci, C2) + s (x, 6, Q, Ci, C2) + {x ^ y) - 8t i-^, /xr 

(41) 

with r (Ci/6, /ir) given by Eq. (|35D and /ir = C4/(x, y)Q. As we shall show below, the effect 
of including the effective quark mass in the hard part of the form factor is almost negligible, 
as one should expect on theoretical grounds. 

Before we go beyond the leading order in the perturbative expansion of the hard- 
scattering amplitude, Th, let us pause for a moment to comment on the pion wave func- 
tion. We have pro-actively indicated in Eq. (|40|) that the asymptotic distribution amplitude 
^as{x) = 6xx will be used. 

A few words about this choice are now in order. 

Hadron wave functions are clearly the essential variables needed to model and describe 
the properties of an intact hadron. In the past, most attempts to improve the theoretical 
predictions for the hard contribution to the pion form factor have consisted of using end- 
point concentrated wave functions (distribution amplitudes). In this analysis we refrain from 
using such distribution amplitudes of the Chernyak-Zhitnitsky (CZ) type referring for a 
compilation of objections and references to |^ (see also [Q), and present instead evidence 
for an alternative source of enhancement due to the non-perturbative power correction in 
the running coupling. 



This IR-enhancement effect was found in to be quite significant, even for the asymp- 
totic solution (to the evolution equation) which has its maximum at x = 1/2. Indeed, 
the IR-enhanced hard contribution can account already at leading perturbative order for 
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a sizable part of the measured magnitude of the electromagnetic pion form factor, though 
agreement with the currently available experimental (low-momentum) data calls for the 
inclusion of the soft, non-factorizing contribution (cf. Eq. (^)) pl|j6^j63|j47[] - even if the 
NLO correction is taken into account (see Sect. ^). Nevertheless, the true pion distribution 
amplitude may well be a "hybrid" of the type •l'^'^"" = 90%$^" + 9%$^^ + l%Cf^^\ where 
the mixing ensures a broader shape with the fourth-order, "Mexican hat" -like, Gegenbauer 
polynomial C4 , being added in order to cancel the dip of $cz at x = 1/2. 

First tasks from instanton-based approaches show that the extracted pion distribution 



amplitudes are very close to, albeit somewhat broader than, the asymptotic form [64,65,66 



Similar results were also obtained using nonlocal condensates The discussion of 

non-asymptotic pion distribution amplitudes will be conducted elsewhere. 



4. Pion form factor to order 



.a" 



Next, we generalize our calculation of the hard contribution to the pion form factor by 
taking into account the perturbative correction to Th of order a^, using the results obtained 
in [^,^,^,^], in combination with our analytical, i.e., IR-finite (IRF) factorization and 
renormalization scheme. 

To be precise, we only include the NLO corrections to Th, leaving NLO corrections to the 
evolution of the pion distribution amplitude aside. The reason is that for the asymptotic 
distribution amplitude, at issue here, these corrections are tiny, appearing first at NLO 
|69| , |27[| . For sub-asymptotic distribution amplitudes, however, evolutional corrections ||69 



have to be taken into account. Strictly speaking, the calculation below is incomplete, the 
reason being that the transverse degrees of freedom in the NLO terms of Th have been 
neglected, albeit the intrinsic ones in the wave functions have been taken into account - in 
contrast to other approaches [^]. Hence, our prediction should be regarded rather as an 
upper limit for the size of the hard contribution to the pion form factor than as an exact 
result. Taking into account the fc_L-dependence of Th at NLO, as we did for the leading part, 
this result might be somewhat reduced as shown for the pion in I^O] and for the nucleon in 
ITOI (for a comprehensive discussion of k± effects, we refer to [0), though we expect that 
due to IR-finiteness, this reduction should be rather small and the quality of our predictions 
almost unchanged. Note in this context that we always refer to the asymptotic distribution 
amplitude of the pion. Broadening the pion distribution amplitude would lead to a larger 
normalization of (form-factor) magnitudes. We would also like to emphasize that other 
higher-twist contributions of non-perturbative origin, as those mentioned before, may also 
raise the magnitude of the form factor. However, such contributions are not on the focus of 
the present work. 

Applying these assumptions, Eq. (^0]) extends to NLO to read 



{xx + yy) 



X exp 



exp 



Mxx yy^ 



xK {^/xyQh) exp (-5 (x, ?/, 6, Q, Ci, C2, C4)) 
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1 + ^ (/uv (a;,2/,QV/^R) + /iR (a;,i/,<5V/^F) + /c(a;,i/) 



(42) 



where the Sudakov form factor, including evolution, is given by Eq. (^), = Ci/b, and 
the functions fi are taken from They are given by 



/uv {x,y,Q^/fJ^R 



^ (3 + In (xy)) Q In (xy) - In ^ 

— [—34 + 12 In {xy) + In x In ?/ 

+ In X In ?/ — In X In 1/ — In x In y 
+ (1 - X - ?/)i^(x, y) + R{x, y)] 



(43) 



and are related to UV and IR poles, as indicated by corresponding subscripts, that have 
been removed by dimensional regularization along with the associated constants ln(47r) — 7e, 
whereas fc{x,y) is scale-independent. In evaluating expression fc in (|43|), we found it 
particularly convenient to use the representation of the function H{x,y) given by Braaten 
and Tse |g. 



H{x,y) 



1 — X — y 



[Li2 (x) + Li2 {y) - Li2 (x) - Li2 (y ) + In X In ?/ - In X In y] 



(44) 



where again Li2 denotes the dilogarithm function. For the function i?(x, y) we have used 
the expression derived by Field et al. [^], except at point x ^ y, where we employed the 
Taylor expansion displayed below: 



R{x,y) 



3{-l+y) y^ 



x-y 



3(_l+^)2^3 
'30(-l + y)3^4 



"(-1 + 33?/ -45 + 13 2/3) Iny 
+ y {-l+y + {9-13y)y\ny) 
j-l+2/)2 (-1 + 162/) Iny 
+ y (^-l + 13y ~12y^ + 2y^ Iny 

{-l + yf (9- 148?/ + 9 2/2) In^ 
- 2/ (9 - 148 2/ + 328 2/2 - 189 y^ + y^{5 + 9 y) In y 



(45) 



Note that this expression does not reproduce its counterpart in |Q. It must be remarked 
once again that evaluating Eq. (|i2| ) there is an analytization ambiguity similar to that 
encountered in the calculation of the Sudakov exponent. This question will be addressed 
elsewhere. 

Having developed in detail the theoretical apparatus, let us now turn to the concrete 
(numerical) calculation of the pion form factor at NLO. 



21 



5. Numerical analysis 



This section implements factorization on the numerical level, thus providing the bridge 
between the analytic framework, developed and discussed in the previous sections, and 
numerical calculations to follow in the subsequent section. This is done by appropriately 
defining the accessible phase space regions (kinematic integrals) of transverse momenta (or 
equivalently transverse distances b) for gluon emission in each regime, making explicit the 
inherent kinematical restrictions on the momenta of hard (soft) gluons due to factorization. 
The numerical analysis below updates and generalizes our previous investigations in Refs. 



In order to set up a reliable algorithm for the numerical evaluation of the expressions 
presented above, we have to ensure that this is done in kinematic regions where use of fixed- 
order or re-summed perturbation theory is legal. Further, expedient restrictions have to 
be imposed to avoid double counting of gluon corrections by carefully defining the validity 
domain of each contribution to the pion form factor, in correspondence with Fig. |^. These 
kinematic constraints are compiled below. 

Kinematic cuts 

1. Ci/b > Aqcd; otherwise the whole Sudakov exponent exp(— S") (cf. Eq. (|4T|) ) is con- 
tinued to zero because this large-6 region is properly taken into account in the wave 
functions. This condition excludes from the re-summed perturbation theory soft glu- 
ons with wavelengths larger than Ci/A, which should be treated non-perturbatively. 
In other words, it ensures the separation (factorization) of the effective sub-sector from 
the genuine non-perturbative regime (cf. Fig. |^). 



2. C2^Q > Ci/b] otherwise each Sudakov exponent exp [— s (^, 6, Q, Ci, C2)] in Eq. (^T]) 



is "frozen" to unity because this small-6 region is dominated by low orders of per- 
turbation theory rather than by the re-summed perturbation series, and consequently 
contributions in this region should be ascribed to higher-order corrections to Th, which 
we have taken into account explicitly at NLO. This condition establishes proper factor- 
ization between re-summed and fixed-order perturbation theory and helps avoid double 
counting of such contributions (always working in the gauge y4+ = 0). Yet evolution 
is taken into account to match the scales in our "gliding" factorization scheme. 



3. C4/(x, y)Q > Ci/b; otherwise the evolution time r {Ci/b, /ir) in Eq. ( PD is contracted 
to zero, i.e., evolution is "frozen". The renormalization scale should be at least equal 
to the factorization scale, so that the running coupling has always arguments in the 
range controlled by (re-summed or fixed-order) perturbation theory. 

4. C4f{x,y)Q > C2iQ', otherwise evolution to that scale is "frozen" because this region 
is appropriately accounted for by the Sudakov contribution. This helps avoiding dou- 
ble counting of terms which belong to the re-summed rather than to the fixed-order 
perturbation theory. (No overlap at the boundary characterized by the scale /i^ in 
Fig. I). 
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5. Cif{x,y)Q > Ci/b; otherwise the two scales /iR = /ip = Ci/b are identified in the 
function fuY{x,y) (by the same reasoning as above). If /ir < Aqcd, then f\j\{x,y) is 
set equal to zero. 

6. Ci/h > Aqcd; otherwise the function fin{x,y) is set equal to zero. The last two 
restrictions exclude contributions from perturbative terms when they are evaluated in 
the non-perturbative kinematic domain. 



To illustrate the difference in technology between approaches employing the conventional 
expression for the full Sudakov exponent [^,^, on one hand, and our analysis, on the 
other, we show exp(— S") graphically in Fig. |^ for three different values of the momentum 



transfer and ^ = 1/2. In contrast to Li and Sterman [^, the evolutional contribution 
is not cut-off at unity, whenever C2^Q < Ci/b. The dotted curve shows the result for 
Eq. (|39D without this cutoff. One infers from this figure that their suggestion to ignore 
the enhancement due to the anomalous dimension does not apply in our case because the 
IR-modified Sudakov form factor is not so rapidly decreasing as b increases, owing to the 
IR-finiteness of a^'^. Indeed, as Q becomes smaller, exp(— S") remains constant and fixed to 
unity for increasing b, providing enhancement only in the large-6 region before it reaches 
the kinematic boundary Ci/b = Aqcd, where it is set equal to zero. As a result, for 
small Q-values, like Qi = 2 GeV, the enhancement due to the quark anomalous dimension 
cannot be associated with higher-order corrections to Th, since it operates at larger b- 
values, and for that reason it should be taken into account in the Sudakov contribution. 
Only for asymptotically large Q values, when the IR-modified Sudakov form factor and the 
conventional one become indistinguishable, the evolutional enhancement becomes a small 
effect - strictly confined in the small-6 region - and can be safely ignored. On the other 
hand, because the Sudakov exponent is bounded at fixed Q^, the Sudakov exponential 
remains finite until the edge of phase space, Ci/ba < Aqcd, also providing IR enhancement. 
This behavior is best appreciated by comparing the dashed and dotted curves, both at 
= 5 GeV, in Fig. |. 

The behavior of the Sudakov form factor, we stress, shows that power-induced sub-leading 
logarithmic corrections are relevant in the range of currently probed momentum-transfer 
values. Hence, the advantage of employing such a scheme to calculate hadronic observables, 
for instance the pion form factor, is that the hard (perturbative) contribution is enhanced, 
relative to the calculation in the MS scheme, and the self-consistency of the perturbative 
treatment towards lower values, where it is not justified, is significantly improved (no 
enhancement caused by the Landau obstruction; hence better scaling). 

This is because the range in which soft gluons build up the Sudakov form factor is 
enlarged and inhibition of bremsstrahlung sets in at larger Q^. Let us mention in this context 
that power contributions in the radiative corrections to the meson wave function could lead to 
suppression of soft gluon emission at large transverse distance b. Indeed, Akhoury, Sincovics, 
and Sotiropoulos |]72| have re-summed such power corrections, associated to IR renormalons, 
with the aid of an effective gluon mass. They found Sudakov-type suppression on top of 
the Sudakov suppression discussed so far. The discussion of such IR-renormalon-based 
contributions in conjunction with our IR-finite approach will be presented elsewhere. 
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Figure 3. Behavior of the Sudakov form factor with respect to the transverse separation h for three 
representative values of the momentum transfer Q^: Qi = 2 GeV, Q2 = 5GeV, and Q3 = 10 GeV, 
with all = 1/2, and where we have set Ci = 2e~'^^, C2 = 6"^/^ and Aqcd = 0.242 GeV. The 
dotted curve shows the result obtained with a^^, and Aqcd = 0.2 GeV for Q2 = 5 GeV, using the 
same set of Ci as before. Notice that in this case, evolution is limited by the (renormalization) scale 
/iR = t = {maxy^xy Q, Ci/6}, as proposed in |29|. However, the enhancement at small 6- values due to 
the quark anomalous dimension is not neglected here. 
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Table 2. Different sets of coefficients Cj and values of the i^-factor and the quantity 
(cf. Eq. (^)) K = In (C^ / 4lC2) , corresponding to different factorization and renormalization 

prescriptions. The choice of C4 in this work corresponds to BLM-type commensurate scale fixing. 



Choice Ci 


Scheme pai 

c. = j,c^'' 1 


•ar 
42 


neters Q 
1 Cs 


C4 K 




canonical 2 exp ( — 7e ) 
SSK |22] exp -i (27E - 1) 
this work 2 exp (— 7e) 


I/V2 
exp (-1/2) 




2 exp (-7e) 
exp -i (27E - 1) 
2exp(-7E) 


4.565 
2.827 

exp (-4/3) 4.565 


-0.307 





6. Phenomenology 



Let us now present phenomenological applications of our scheme. Using the techniques 
discussed above, we obtain for the electromagnetic pion form factor the theoretical predic- 
tions shown in Fig. ^. A set of constants Ci, {i = 1, 2, 3) which eliminate artifacts of dimen- 
sional regularization, while practically preserving the matching between the re-summed and 
the fixed-order calculation, are given in Table |] in comparison with other common choices 
of these constants. Moreover, this factorization scale setting enables us to naturally link 



our scheme to the BLM commensurate scale method in fixing the renormalization scale. 
Indeed, since the adopted value of Ci eliminates both the log term in the ii'-factor (see 
Eq. (|3^), which contains the /^-function, and also the scheme-dependent term B in the 
cusp anomalous dimension (see Eq. (|5BD), this choice corresponds to a conformally invari- 
ant framework with /5o = 0, and therefore connects to the commensurate scale procedure. 
Hence, we set C4 = C2exp (—5/6), which, for our choice of C2 = exp (—1/2), rescales Q in 
the MS scheme, we use, by a factor of exp (—4/3). In addition, to avoid large kinematical 
corrections due to soft gluon emission, we set f{x,y) = y/xy to link the renormalization 



scale to the typical momentum flow in the gluon propagators ||5^. In this way, scheme and 
renormalization scale ambiguities are considerably reduced, as the theoretical predictions 
are evaluated at a physical momentum scale: 

AiBLM = /iRexp(-5/6) , (46) 

where 



/iR = CJ{x, y)Q = C^y/^Q . (47) 

We emphasize, however, that these favored values of the scheme constants by no means 
restrict the validity of our numerical analysis. They merely indicate the anticipated appro- 
priate choice of the factorization and renormalization scales with respect to observables and 
theoretical self-consistency. Other choices of these parameters do not change the qualitative 
features of our predictions. 

Before we proceed with the discussion of these results, let us first present the theoretical 
prediction for the pion-photon transition form factor F^^*^(Q^,g^ = 0) in which one of the 
photons is highly off-shell and the other one is close to its mass-shell. In leading perturbative 
order, this form factor is given by the expression (cf. [0) 
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Figure 4. Space-like pion form factor calculated within our IR-finite scheme with <I>as and including 
an effective (constituent-like) quark mass of rriq = 0.33 GeV in the pion wave function. The broken lines 
show the IR-enhanced hard contributions obtained with our scheme using commensurate scale setting; LO 
calculation (dashed-dotted line); NLO calculation (dotted line). The dashed line gives the result for the 
soft, Feynman-type contribution, computed with m-q = 0.33 GeV in the pion wave function. The solid line 
represents the sum of the NLO hard contribution and the soft one. The data are taken from [[73l|74|] . 



26 



0.3 



CLEO 
• CELLO 




0.0 2.5 5.0 7.5 10.0 

Q' [GeV^] 

Figure 5. Pion-photon transition form factor for the asymptotic distribution amplitude, calculated 
in the IRF scheme. The dashed line shows the prediction for the BHL ansatz for the pion wave func- 
tion which includes an effective (constituent-like) quark mass of m-q = 0.33 GeV. Commensurate scale 
setting is used, i.e., C4 = C2exp(— 5/6) with C2 = exp(— 1/2). For comparison, we also show the 
prediction (solid line) obtained without a quark mass and using a non-commensurate renormalization scale 
(C4 = C2 = exp(— 1/2)). The horizontal line represents the asymptotic behavior. The data are taken 
from (7|@. 

F., (q^) = -|- I' dx f dhh^J^l^^ exp {-x-xh^/^^l) 

X exp (-/^as^^) (4vri^o (Vi hQ)) e'^^^ , (48) 
where the Sudakov exponent, including evolution, has the form 

S^^{x,x,h,Q,Ci,C2,C4) = s{x,h,Q) + s{x,h,Q) - At {^,11^ . (49) 

The main difference relative to the previous case is that this form factor contains only one 
pion wave function, whereas the associated hard-scattering part, being purely electromag- 
netic at this order, does not depend directly on Og. The only dependence on the (running) 
strong coupling enters through the anomalous dimensions in the Sudakov form factor. The 
result of this calculation is displayed in Fig. |^. 

All constraints on kinematics set forward in the numerical evaluation of the electro- 
magnetic pion form factor are relevant to this case too, except the requirement which 
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deals specifically with the choice of the renormalization scale, which now is set equal to 
//r = CixQ because only one pion wave function is involved. Another reasonable choice 
would be /iR = C^\fxi,Q^ which entails evolution to a lower scale, hence reducing evolu- 
tional enhancement through r /zr) by approximately 6%. 
Let us now discuss these effects more systematically. 

It is obvious from Fig. | that the IR-enhanced hard contribution to F^(Q^) with opti- 
mized choice of scales is providing a sizeable fraction of the magnitude of the form factor 
- especially at NLO. This behavior is IR stable from low to high values, exhibiting al- 
most exact scaling (in accordance with the nominal scaling of the leading twist prediction), 
which shows that the analytic coupling is almost constant in a wide range of values. 
In contrast to other approaches [|3ll , p7i^ , which involve a running coupling without an 
IR-fixed point, there is no artificial rising at low of the of the hard form factor, resulting 
from the unphysical Landau pole. Furthermore, by employing a commensurate scale setting 
to fix the renormalization point, the scheme and renormalization-prescription dependence 
of our predictions has been minimized. In addition, the imposed kinematical constraints in 
our numerical analysis ensure that the contributions, originating from different phase space 
regions, do not overlap to give rise to double counting. 

The reduced sensitivity of the perturbatively calculated hard form factor to the endpoint 
region, is also reflected in its saturation behavior. One sees from Fig. ^ that the bulk of 
the scaled form factor Q^F^ (Q^) is already accumulated below fecrAqcD/C'i < 0.5, i.e., for 
short transverse distances, where the application of perturbative QCD is self-consistent. All 
curves shown rise very steeply to their full height at the integration cutoff 6cr = Ci/Aqcd, 
beyond which they flatten out, indicating that remaining contributions are truly of nonper- 
turbative origin. One observes that the perturbative treatment starts to be reliable already 
at = 4 GeV^ and improves further, albeit not dramatically, as the momentum transfer 
increases. A fast saturation behavior in the small 6-region, where contamination with non- 
perturbative contributions is still not serious and the coupling constant is small, is considered 
as a standard to judge the self-consistency of the perturbative method applied. Though the 
Sudakov form factor contains considerable contributions from gluons with wave lengths of 
the order of Ci/Aqcd due to the IR finiteness of the running coupling (see Fig. ^ - espe- 
cially at low momentum transfer - we realize that the form factor itself, i.e., the physical 
observable, does not receive strong contributions from this endpoint (6) region. Moreover, 
the form factor calculation does not receive large contributions from the endpoint region in 
X as well, as we use only the asymptotic pion distribution amplitude (or such close to it). 
Hence, from the theoretical point of view, the quality and self-consistency of the perturbative 
treatment have been improved relative to previous approaches PU| , PU |. 

Figure |^ shows the influence of the effective quark mass on the pion form factor. The 
designations are as follows: The solid line plots (P'^lijF'^l^^ f^^^ and the dotted hne 

[f^^^^qI F^^^^^\^^^^ . It is clearly obvious that the effect of the quark mass on the hard 
part is negligible in size and does not depend on the variation in Q^, whereas the soft contri- 
bution gets significantly reduced as grows. The dashed line, standing for the expression 
(-^comm/-^non-^ Same figurc quantifies the effect of using a commensurate 

scale setting for the renormalization scale. As one sees, this amounts to an enhancement 
factor of about 1.5. 

The advantages of our framework may become more transparent by comparing our results 
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Figure 6. Saturation behavior of the pion electromagnetic form factor, calculated in the IRF scheme 
at NLO with commensurate scale setting and including a mass term (with mq = 0.33 GeV) in the BHL 
ansatz for the soft pion wave function. The scheme parameters are defined in Table |[ Here bcr denotes the 
integration cutoff over transverse distances in Eq. E3. The momentum transfer values are as in Fig. H. 



Table 3. Calculated pion form factor at two values of Q^. The first two columns show the 
results obtained in the present work in comparison with those calculated by Jakob and Kroll (JK) 
|l30|| (third column), and by Melic, Nizic, and Passek (MNP) ||2^ (last two columns). 



g2 [GeV2] 


LO 


LO+NLO 


LO 


LO 


LO+NLO 




(this work) 


(this work) 


(JK) 


(MNP) 


(MNP) 


4 


0.128 


0.191 


0.08 


0.131 


0.211 


10 


0.137 


0.186 


0.08 


0.109 


0.164 



with those obtained in other analyses. This is done for the pion form factor in Table ^ 



Comparison of our values with those calculated by Jakob and Kroll shows that the 
suppression of the hard part of the form factor due to the inclusion of transverse degrees 
of freedom is counteracted by the power-induced enhancement, amounting to an average 
enhancement of about 50% relative to their values. This is achieved by using (almost) 

the same root mean square transverse momentum of (k^^)^ ^ = 0.352 GeV, as in their 
analysis, and with a reasonable probability for the valence Fock state of Pqq = 0.306 (see 
Table |I]). The inclusion of an effective (constituent-like) quark mass in the Gaussian ansatz 
for the distribution of intrinsic transverse momentum in the pion wave function changes 
dramatically the fall-off behavior of the soft contribution to the form factor, as compared 
to the JK analysis, though its maximum size remains almost unchanged, and its influence 
on the hard part is very small (cf. Figs. ^ and |^). Indeed, one infers from Fig. | that 
becomes equal to F^^'''^ already at Q'^ ~ 18 GeV^ (LO result), or even at ~ 12 GeV^ 
when the NLO corrections are included. This behavior of falls well in line with the 
correct behavior of vl>soft ^^j. = q g^nd —oo, restored by the mass term and the 
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Figure 7. Effect of using a BHL-type ansatz for the pion wave function with an effective (con- 
stituent-like) quark mass of mq = 0.33 GeV. The sohd hne plots the ratio (-?*'mq=o/-^mq'^o) Iramm ^^'^ 



the dotted line the ratio (^F^ 

the ratio ( -^^commZ-^non— comm 



soft 

'q = 



'Soft 



iNLO 



mqT^o) I 



versus the momentum transfer Q . The dashed line shows 



of the hard part of the pion form factor for a BLM commensurate 



scale setting relative to a conventional one with C4 = C2 = exp (—1/2) (cf. Table g). 
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arguments on the nonperturbative vacuum dynamics given above. 

On the other hand, comparison with the values computed by Mehc, Nizic, and Passek 
p7| at leading order, by completely ignoring transverse degrees of freedom, reveals that 
in the domain, where the influence of the Landau singularity has died out (values for 

= 10 GeV^ in Table there is still enhancement of about 17%. Comparing our results 
with theirs at next-to-leading order, we conclude that our choice of scheme and renormal- 
ization scales is consistent with a proper matching between gluon corrections, calculated 
on a term-by-term perturbation expansion (NLO corrections to Tfj), and those due to the 
re-summed perturbative series (Sudakov form factor). Therefore, double counting of such 
contributions in our scheme, if any, must indeed be negligible. Moreover, the scaling behav- 
ior of the calculated perturbative (hard) form factor is considerably improved, complying 
with the nominal scaling of the leading twist prediction. Indeed, one observes (cf. ^ that the 
deviation from exact scaling, associated with NLO evolutional corrections of the asymptotic 
distribution amplitude, is, as stated before, negligible. 

The illustration of the enhanced form-factor behavior (always assuming the asymptotic 
form of the pion distribution amplitude) is given in Fig. |^ in terms of the ratio between 
F^SK (Q2^^ calculated in this work, and F^^^ {Q"^), obtained by Mehc et al. in One 
sees from that figure that at values up to about 5 GeV^, this ratio is less than unity, clearly 
exhibiting the singular IR-behavior of the conventional as(Q^) representation, employed by 
these authors. Contrary to that, above approximately 10 GeV^, this ratio scales with at 
a fixed value of about 1.25. Hence, restoring analyticity of the effective QCD coupling (by 
a power correction term), removes the artificial raise of the form factor, owing to the rapid 
increase of the perturbative coupling at low momentum, and stabilizes its low-Q^ behavior, 
providing enhancement only in that momentum region which is controlled by self-consistent 
perturbation theory. 

Let us turn again to the calculation of F-^^^Q"^). Figure || shows our theoretical pre- 
dictions for this form factor using the same set of scheme parameters Ci,C2,C3, given in 
Table The dashed line includes a quark mass term and employs commensurate scale 
setting for the renormalization point. The solid line shows the prediction for mq = and 
a non-commensurate renormalization scale, with (74 = 6*2 = exp(— 1/2). This latter curve 



reproduces the recent high-precision CLEO ||75[ and also the earlier CELLO |^ data with 
almost the same numerical accuracy as the dipole interpolation formula. However, we regard 
the lower curve as being more realistic because a physical renormalization scale has been 
used (provided our choice of C4 = C2exp(— 1/2) is approximately correct). Remarkably, 
the predicted magnitude of Q^F^^, being somewhat below the data, allows some broadening 
of the pion distribution amplitude, as recently found in instanton-based approaches |51,5B[| 
or using non-local condensates |]67| , |68| . 

The sensitivity of the pion-photon transition form factor to the quark mass and the 
commensurate scale setting is discussed in Fig. Rl. 



Our prediction is consistent with the result obtained by Brodsky et al. in IT^l, who also 
use commensurate scale setting and include in addition the LO QCD radiative correction 
to Ftt^ with a running coupling "frozen" at low momenta by virtue of an effective gluon 
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Figure 8. The ratio R-^ (Q^) = F^SK (gS-j /pMNP (-g2-) ^^^^^^ momentum transfer Q^. F^^^ is 
the pion form factor given in Eq. dH), F^NP ^j^g expression derived in |7| (cf. Table |), using the asymp- 
totic pion distribution ampUtude. The decrease of this ratio below ~ 5 GeV'^ signals the breakdown 
of perturbation theory in the calculation of |^ owing to the Landau singularity in the conventional ttg 
representation they use. 
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[GeV^] 



Figure 9. Effects of commensurate scale fixing and a non-zero quark mass in the BHL ansatz 
for the pion wave function in the pion-photon transition form factor. The solid line shows the ra- 
tio ^F^"^ '^7-^^'^^'^) Icomm ^'^^ ^ commensurate scale fixing and the dotted line the same ratio for 

a non-commcnsuratc scale fixing: (^F^ °/-^^'''^°) Inon-comm with C4 = C2 = exp(— 1/2). The 
dashed line effects the difference between using a commensurate scale fixing and a conventional one with 
C4 = C2 = exp(-l/2) for the ratio L^^O with a non- vanishing quark mass 

m„ = 0.33 GeV. 
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t igure 10. The ratio R^i (Q^) = (q2) /qS^ (q2) 

versus the momentum transfer Q . is 
the analytic running couphng in one- loop approximation (see, Eq. (|l|)). The effective charge cty, used by 
Brodsky et al. in ||7^, is defined in Eq. (jsj). For more explanations and the definition of i?„2 (Q'^), see in 
the text. 
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Table 4. Values of the scaled space-like pion form factor, calculated in our IRF scheme at 
different momentum transfers . Q^F^^ (Q^) is the LO result given by Eq. ( |40| ) and represented 
by the dashed-dotted line in Fig. |. Q^^nlo ^q2^ (dotted hne in Fig. ||) is the expression displayed 
in Eq. ( |4^ ) and comprises the LO and NLO contributions to the hard-scattering part (for more 
details, see Sect. Q). These results were obtained with a non- factorizing BHL-type ansatz for the 
pion wave function (i.e., with an effective (constituent-like) quark mass mq = 0.33 GeV in the pion 
wave function), and employing BLM commensurate (renormalization) scale setting. The last two 
columns show the results for the pion-photon transition. Q^F^^ {Q^) stands for the expression (48) 
and commensurate scale setting, whereas Q'^F^~^ (Q^) shows the results without the inclusion 
of a quark mass and with a non-commensurate renormalization scale (cf. dashed and solid lines in 
Fig. |5|, respectively). The asymptotic pion distribution amplitude is always assumed. 



Q2 [GeV2] 


Q'F^^ (Q2) 


g2^NLO (q2) 


Q'F^, (Q2) 




2 


0.1121 


0.1831 


0.1180 


0.1370 


4 


0.1282 


0.1907 


0.1317 


0.1576 


6 


0.1340 


0.1904 


0.1375 


0.1668 


8 


0.1364 


0.1882 


0.1407 


0.1721 


10 


0.1373 


0.1856 


0.1428 


0.1755 


15 


0.1368 


0.1790 


0.1458 


0.1803 


20 


0.1351 


0.1731 


0.1474 


0.1828 


30 


0.1312 


0.1639 


0.1491 


0.1854 


40 


0.1275 


0.1568 


0.1500 


0.1867 



mass.u The close resemblance between the two approaches becomes apparent by comparing 
the corresponding running couplings against the momentum transfer. In Figure |10| we show 
the ratio (solid line) 



with {Q^) given by Eq. (P and the coupling (effective charge) in the so-called V scheme, 
defined by 



where Ay = 0.16 fT^] and = 0.19 GeV^ The dashed line {Ra,2 (Q^)) represents this 
ratio with Ay set equal to A = 0.242 GeV in Eq. ([51|). Though, strictly speaking, it is 
inconsistent to equalize scheme-dependent parameters, the message of this figure is that the 
two parameterizations are very close to each other, albeit the analytic coupling has a larger 
normalization at low Q^. 



The connection between the modified convolution scheme, which explicitly retains transverse 
degrees of freedom, and the use of an effective gluon mass to simulate the effect of the Sudakov 
suppression factor, was discussed in |7l| . 
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Closing our discussion of the photon to pion transition, let us mention that other authors 
7^,79.8^ obtain similarly good numerical agreement of Q^F^^* (Q^) with the experimental 



data, following different premises based on QCD sum rules. 

Finally to facilitate a more detailed comparison of our results with other approaches and 
experimental data, we compile in Table ^ the obtained values of the (scaled) pion electromag- 
netic (LO and NLO) and photon to pion (LO) form factors at different momentum transfers 
. These form factors are calculated with a non-factorizing BHL-type ansatz for the pion 
wave function (hence including an effective quark mass in the Gaussian distribution for the 
intrinsic transverse momentum) and using a BLM commensurate fixing of the renormaliza- 
tion scale. In the case of the pion-photon transition form factor, we also show the result 
setting the constituent quark mass equal to zero and employing a non-commensurate renor- 
malization scale - in analogy to our previous analysis in |]22|. In all cases, the asymptotic 
form of the pion distribution amplitude is assumed. 



7. Summary and conclusions 

Let us summarize the hallmarks of the presented methodology. We have developed in 
detail a theoretical framework which self-consistently incorporates effects resulting from a 
modification of the strong running coupling by a non-perturbative minimum power cor- 
rection IQ which provides IR universality. Though a deep physical understanding of such 
contributions is still lacking, we have given, as a matter of practice, quantitative evidence 
that using such an analytic running coupling it is possible to get an IR-enhanced hard con- 
tribution to the electromagnetic form factor F-j^lQ"^) by employing only asymptotic (like) 
forms of the pion distribution amplitude, hence without recourse to end-point concentrated 
distribution amplitudes. 

The presented IR-finite factorization and renormalization scheme makes it possible to 
take into account transverse degrees of freedom both in the pion wave function as well 
as in the form of Sudakov damping factors , without entailing suppression of the (pion) 
form-factor magnitude resulting from severe IR regularization. In addition, use of this mod- 
ified form of as(Q^) renders the theoretical predictions insensitive to its variation with at 
small momentum values, thus remarkably improving their scaling behavior, in accordance to 
the nominal scaling of the leading twist contribution. Similarly, the saturation behavior of 
the pion form factor (versus the impact separation) is significantly improved and the scaled 
hard form factor reaches much faster a plateau, accumulating its magnitude in the region 
of small transverse distances where use of perturbation theory is legal. An appropriate 
choice of the factorization (scheme) scales and the strict separation between gluonic contri- 
butions from fixed-order and re-summed perturbation theory helps avoid double counting 
of higher-order corrections, enforcing this way the self-consistency of the whole perturbative 
treatment in a wide range of momentum transfer. Moreover, adopting the BLM commensu- 
rate procedure in order to choose an optimized renormalization scale, and thus minimize the 
renormalization scheme dependence, we have calculated the pion form factor including the 
NLO radiative correction to the hard-scattering amplitude. In contrast to other approaches, 
we employ a BHL-type of ansatz for the distribution of the intrinsic transverse momentum 
in the pion wave function which includes a mass term. This term, resulting from the non- 
perturbative QCD vacuum structure, ensures suppression of (x, k_|_) for k_|_ = and 
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ks — > — oo. and yields to a stronger fall-off of the soft, non-factorizing contribution to the 
form factor at momentum-transfer values around 20 GeV^ and beyond. Hence, the leading- 
twist predictions of QCD are remarkably confirmed at still higher Q^, whereas at lower 
momentum values Feynman-type contributions dominate. In this region other higher-twist 
contributions may also be important. 

The same procedure applied, without any scheme parameter re-tuning, to Fj^o^*^ yields 
a prediction which is consistent with, though somehow below, the experimental data of the 
CLEO and CELLO collaborations, and allows therefore for a mild broadening of the (true) 
pion wave function, as indicated by instanton-based approaches. 

We believe that the insight gained through our analysis gives a strong argument that a 
power correction in the running coupling of QCD, as proposed by Shirkov and Solovtsov, has 
important consequences and provides a convenient tool to improve theoretical predictions 
based on perturbation theory. 
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